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Abstract 

The detection and characterization of a non-classical-squeezed state of light, by using two dif¬ 
ferent schemes, will be presented . In the first one, in an one-dimensional cavity with moving 
mirror (non-stationary Casimir effect) in the principal mode, we study the photon creation rate 
for two modes (squeezed and coherent state) of driving field. Since the cavity with moving mirror 
(similar to an optomechanical system) can be considered analogue to a Kerr-like medium, so that 
in the second scheme, the probability amplitude for multi-photon absorption in a nonlinear (Kerr) 
medium will be quantum mechanically calculated. It is shown that because of presence of nonlinear 
effects, the responses of these two systems to the squeezed versus coherent state are considerably 
distinguishable. The drastic difference between the results of these two states of light can be viewed 
as a proposal for detecting of non-classical states. 

PACS numbers: 42.50, 42.65 


* r.roknizadeh@gmail.com 


1 



I. INTRODUCTION 


In recent years many experimental and theoretical progress in high-Q electromagnetic 
cavities have been made. Among them, there is a phenomenon that can be seen in the name 
of Non-Stationary Casimir Effect (NSCE) or Dynamical Casimir Effect (DCE) caused by 
changes of (vacuum) quantum states of holds due to a fast time variations of positions (or 
properties) of boundaries confining the fields [1-3]. The DCE was invented by Yablonovitch 
[4] and Schwinger [5] for the first time. The DCE has been observed in many systems 
such as: single mirror cavity [6], cavities with moving mirrors [1-3, 7], circuits [8] and 
parametrically excited media [9]. Amplification of quantum (vacuum) fluctuations in these 
macroscopic systems is a common feature. In cavity with moving mirrors, DCE concerns 
to the reconstruction of the quantum state of a held, due to a time dependent geometrical 
configuration. Many phenomena [10] such as photon creation [7], interaction between the 
held and detectors [7, 11, 12], generation of nonclassical (in particular squeezed) states [13], 
etc. are observed in these cavities. Dodonov et al have investigated the photon creation in 
one-dimensional cavity in vacuum state (empty cavity) [7]. 

Sousa et al [14] recently have shown that the nonlinear Kerr effect can be observed in 
cavities with DCE. In other words, a realistic description of cavity DCE include the Kerr 
nonlinearity. 

In other side, optomechanics as a new held has some common features with DCE (which 
has been observed in cavity with moving mirror). Gong et al , recently have showed that, 
the Kerr nonlinear optical effect can be observed in optomechanical system in vacuum. So 
the cavity with moving mirror can be supposed as a Kerr-like nonlinear medium [15]. 

There are many standard methods for generation and detection of non-classical, in particular 
squeezed, light. The squeezed light detection has been done in homodyne systems [16-19], 
which is a phase-sensitive measurement of the optical held quadratures. The squeezed light 
interaction with atomic systems has been predicted in Refs. [20, 21], While the most 
important advantage of cavity with moving wall (optomechanical system) is detection the 
quantum mechanical effects (such as squeezing) macroscopically. 

Recently, some systems such as non-linear media [22-24] and optomechanical systems 
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[25, 26] and resonators with moving wall [13] have been considered for generating squeezed 
light. But, in nonlinear media with Optical Parametric Oscillator (OPO), the squeezed 
parameters are affected by various experimental features [27-29]. 

Moreover nonlinear processes such as multiphoton absorption can change the photon 
statistics of the incident light and generate a narrower photon distribution than poisson 
distribution [30]. On the other hand, multiphoton absorption processes are very sensitive to 
the type of light state and its fluctuations. The influence of nonclassical light fluctuations 
particularly the squeezed states on the two photon absorption has been studied and its 
absorption rate compared with respect to the phase and amplitude of the squeezed beam 

[31] . The two photon interaction of atoms and squeezed light has been considered in Refs. 
[32, 33]. The effect of the absorption process on the squeezed light has been studied in Ref. 

[32] , And the two photons quantum optics has been investigated by [33]. In [34], as well, 
the multiphoton absorption rate is calculated by quantum mechanical perturbation theory 
in a semi-classical approach. 

Wilson et al [8] in their paper has experimentally shown the effect of DCE in supercon¬ 
ducting circuit. They have mentioned that this system operates as a parametric amplifier 
(nonlinear medium) to generate photons, and the mechanism of generating photons is very 
similar to a nonlinear medium. 

In this contribution, as it has been predicted in Ref. [8], we are going to present two 
different, but somehow related, approaches for detecting the nonclassical-squeezed-state of 
light. 

Both approaches concern the nonlinear effects, which are sensitive to different modes of 
incident field. In the first part of the paper the interaction of squeezed state and a ID cavity 
with moving mirror will be investigated. In the second one the interaction of squeezed light 
with nonlinear media is studied. The outcomes of these two schemes will be compared with 
the situation, when the incident field is in a coherent state. The drastic difference between 
two setups (squeezed and coherent states) can be used for characterizing the detection of 
the nonclassical light. 

In the first step, we will study the photon generation rate in the one-dimensional cavity 
with moving wall (in the principal mode), in which the incident field is in coherent or 
squeezed state. 

In the second step, by quantizing the electromagnetic field, we derive probability amplitude 
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for three photon absorption in a nonlinear medium, by using quantum polarization [35, 36]. 

In sec. (II), we review the pioneer work of Dodonov and Klimov in Ref. [7], to summa¬ 
rize the related results for calculating the photon generation rate in a cavity with moving 
boundaries. In sec. (Ill) we derive photon generation rate in the cavity for two modes of 
light (squeezed and coherent light). The results for two different modes are also plotted in 
a few diagrams. The fundamental concepts of interaction of these two modes and nonlinear 
media is explained in sec. (IV). In this section we obtain the rate of three photons absorption 
in a semi-classical approach And in sec.(V), the problem of the last section is studied in full 
quantum mechanical approach. In sec. (VI) we obtain the rate of absorption for two kinds 
of light and compare them. Conclusions are presented in sec. (VII). 


II. REVIEW OF PHOTON GENERATION IN CAVITY WITH MOVING BOUND¬ 
ARY IN VACUUM STATE 


In order the paper be self contained, we review in this section the related results given 
in [7]. We consider a cavity having two infinite ideal plates, one being at rest at x — 0, and 
the other one moving according to a time dependent function L(t) when t > 0. We consider 
an electromagnetic field whose vector potential is directed along z axis. The field operator 
in Heisenberg representation at t ^ 0 (when both the plates were at rest) is A(x,t). 


-j x - 1 nnx ■ 

A in = 2 y —j= sm —— b n exp{—ioj n t) + c.c, 


n 


( 1 ) 


7m 


where b n is the annihilation photon operator and c u n = —. The field operator can be written 


A(x, t) = 2 ^ —(b n ilj n (x, t ) + c.c.), 

v ^ 

r> » 


n = 1,2, 


where 


( 2 ) 


. mix . . . 

yj ( x , t < 0) = sm —— exp{—iuj n t). 

Lo 

We expend the function i/V(a;,f) with respect to instantaneous basis[7] 

Tikx 


^{xA > 0 ) =J2Q^\t)J 


sm 


m L(ty 


( 3 ) 


( 4 ) 
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where, for t > T (T is the time interval that the wall returns to its initial position L$) 

Qk\t) = ( 5 ) 

' and y/)) v are some constant coefficients. At t > T we introduce a new set of physical 
operators a m and a) m . The Bogoliubov transformation relates the two sets of operators (a, A 
and b,w) 


'y ^ {b n cx nm + i > l l f3 nrn ) : 

n 

in which the coefficients are 



We assume that the wall oscillates according to the following form, 


( 6 ) 

(7) 


L(t) = L 0 [ 1 + esin(2cuit)], ui = (8) 

To 

where £ <C 1. We put the solution (6) in wave equation. After some simple algebraic 
manipulations we arrive to the following equations in terms of coefficients £ and ry. 


£iV) = 

<| n+ 2 ) (r) = <r 2 ) (r)-^(r) , n> 3 


(n+2) / \ (n— 2) / c • (n) / x ^ 0 

nr?} (r) = (t) — (r) , n> 6 


where we have introduced new time scale: 


r = -euit. 


(9) 

( 10 ) 

( 11 ) 

( 12 ) 


( 13 ) 


By using the Laplace transformation the solution of the above diffrential equations will be 
achieved, and rj[ l> and £j^ are obtained accordingly, 

e a> = 2 E(k) + kK(k) 

7T 1 + k 
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( 15 ) 


2 E{k) - kK(k) 
vr 1 - k 


where K{k) and E(k ) are complete elliptic integrals of the first and second kind: 


K(k) 



1 

\jl — k 2 sin 2 (cc) 


dct 


(16) 


And 


m 




k 2 sin 2 (a)da 


k = Vl — e 8r , k = \/l — k 2 = e 4r . 


(17) 


(18) 


III. INPUTING THE SQUEEZED AND COHERENT LIGHT 

The amount of photons created in the mth mode equals the average value of the operator 
since just this opretator has a physical meaning at t > T. 


A. Squeezed light 

We assume that the held in the cavity is initially in squeezed coherent state |a, C)& where 
the subscript b indicates this state is defined with respect to the b n . The average of photon 
creation in mode m is given by, 


A m bipt-i CI b^ m dj m |cq C) b- (19) 

By using the Eq. (6), the above relation can be written in terms of b and b' : 

Am = £>(a, C | ^ ^ (Pn'^n'm A ^n'Pn'm) 

n' 

X J2^n a nrn + tinPnm)\ a ,0b 

n 

= b{a, Cl J2^'m a nmblb n + 

n' n 

^nrn^n'rn^n'^n + A'mtwl) |a, Ob- ( 20 ) 
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This expectation value can be calculated by using the squeezing operator according to the 
following relations [40], 


1 a,C) = s (0 1 «), 

(a) 


^( 0^(0 = A - A, 

( 6 ) 


S f (C)fts(C) = A - U * h n , 

(c) 


b n | O') = Q! O'). 

(d) 

( 21 ) 


where s(£) is squeezing operator and /i = cosh|£|,z/ = e l¥, sinh |£|, and £, <p are squeezing 
parameter and its phase respectively. After some calculation Eq. (20) gives rise to, 

N m = 

'^2{\a nm \ 2 [ii 2 \a\ 2 - 2fiRe(u*a 2 ) + |u| 2 |a| 2 + \u \ 2 ] 

n 

+ 2Re(/3 nm a nm )Re[/i 2 a 2 - 2/iv\a\ 2 - /iv + v 2 a* 2 } 

+ \/3nm\ 2 \/J 2 \®\ 2 ~ ^ 2 |tt| 2 + /i 2 — 2/iRe(u* a 2 )}}. (22) 

The photon generation rate in principal cavity mode [with m = 1] is given according to, 

dN 

dr 

n 

— 2/iR e(v*oi 2 ) + |zv| 2 |ct| 2 + |z/| 2 ) 

+ — Re(£i' ) rji l ' > + £[ n V/j" ) )Re( / u 2 a! 2 — 2/j,u\a\ 2 
n 

— fiu + v 2 a*~) + — | (A i2 |°^| 2 — 2yuRe(z/*a 2 ) 

n 

— u 2 \a\ 2 + fi 2 )}. (23) 

where we have used of Eq. (7). The sums in Eq. (23) can be simplified 

oo „ 

E (24) 

77.—1 

oo 0 

E -dAb = -24 1, 4 1) . (25) 

71 

n=l 
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( 26 ) 
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Therefore we obtain the final form of the generation rate, 


dN 

dr 

(—277^£j^)(^ 2 |a| 2 — 2/iRe(V*a 2 ) + M 2 |a| J + \is\' 2 ) 

+ 2Re(—)Re( / u 2 CK 2 — 2jiv\a\ 2 — fiu + v 2 a* 2 ) 
+ (-2^ 1) ^ 1) )(|a| 2 ( / n 2 + \u \ 2 ) - 2jjRe(u*a 2 ) + fi 2 ). 


(27) 


B. Coherent light 


Now we apply a driving field in coherent state, and follow the same procedure as given 
in the above subsection in order to calculate the photon generation rate.The field is consid¬ 
ered to be in the coherent state |a)&, where the subscript b indicates that these states are 
eigenvectors of the annihilation operator b n . The photon number created in the mth mode 
is given by N m , 

-/V m | d) m d m | Q))^, (28) 

and as we have just now noted, 

b | a)b = a | a)b, (29) 

where the expansion in terms of Fock states is well known, 

\ a )b — exp(-^-L) I P)- ( 30 ) 

p=0 * 

The action of a m on this state will be calculated by using the transformation relation in Eq. 

( 6 ), 


dm\o^ b 


exp( 


— a 


^ ^ y ^ i—r {otnm&npy/P \ P 1 ) 


n 


+ Pnm,dnp\/P + 1 I P + 1)) 

= X^ exp (— ^-^-)-T=X a nmVn | n- 1) 
n 2 vn! 

+ V,„ v" + 1 I rc + 1))- 


(31) 








By applying the same procedure we can get the action of a m on |a)b, and accordingly, 


°° ^ ^2 n 

N m = Vexp(- | a | 2 )—j-K m + P„ m {n + 1)), 

n\ 


n =0 


( 32 ) 


and by replacing a nm and /3 nm from Eq. (7), 


n=0 


2n 


N m = X] eX P( _ I “ 1 2 )‘h“( n ^m )2 + )• 


(33) 


The photon generation rate in principal cavity mode [with m — 1] is 


dA 2 x - 

— = exp(— | a | ) 2^ 


00 ntWfW 
2 n z Sl Si 


a 


n =0 
00 


n\ 


I ( \ 12\ In 2 n + 1 *(") .*(") 

+ exp(— | a | ) 2^ a ~Vi Vi 


n =0 


n\ n 


(34) 


After some calculation and using of Eqs. (26)-(29), the Eq. (34) can be rewrittn, 


dN 


| a | z -2 
+ 2Re(-ryi 1)2 -«e( 1)2 )Re(a 2 ). 


(i)t(i) 


(35) 


Now, we can compare the photon generation rate for the cavity mode in coherent and 
squeezed state, according to the Eq. (27) and Eq. (35), where the explicit forms of rj[ l> and 
as given in relations (14)-(18), are to be used. It should be considered that the intensity 
of coherent and squeezed light is the same and equals to 7 (The coherent state amplitude 
,a, is a complex number). 

Their plots are shown in Fig. 1. 

The drastic difference between different squeezed lights and coherent light in photon 
generation rate may be easily seen. The effect of squeezing parameter has been shown in the 
same figure. The photon generation rate increases with increasing of squeezing parameter. 
This point is considerable that as usual, the squeezing phenomenon occurs in special region, 
so we restrict ourself in the related region [40]. 

It is required to investigate Eq. (27) in some limits. I 11 the limit r —y 0 the variables 
k — » 0 , k —> 1, {E(k),K(k)} —> r/J 1 ’ —» 0 , and —» 1 . The photon generation 

rate in the limit r —» 0 for squeezed state (a 7 ^ 0 ) has the form, 

= —2Re(/i 2 a 2 — 2/iz/|a| 2 — /i/y + v 2 a* 2 ) for r —» 0, (36) 
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FIG. 1. Photon generation rate versus the dimensionless time t for squeezed light with various 
squeezed parameters (|£| = 0.8, |£| = 1, |C| = 1-2 and |C| = 1-5 , = 0.) and coherent light 

(l «| 2 = 7) [26], 



FIG. 2. Plot of -| r= o as a function of the squeezing parameter |G for different values of real 

ar 

coherent state amplitude a = 0 . 2 , 0 . 8 , 2 . 0 , \/7. 


which , for a fixed real value of a, is a monotonic function of |£| (Fig 2). 

The photon generation rate (for squeezed and coherent states) for real value of a is plotted 

in Fig. 3 which concides with Eq. (36) and Fig. 2. 

For completeness the the amount of the photon creation, when the driving field are in 

vacuum | 0 )& and squeezed vacuum states( with different squeezing parameters), is plotted 

in Fig. 4. The relation for vacuum state in this case coincides with the Eq. (6.1) in [7]. 

According to the Eq. (36), —— | T= o is equal to 2/i/y which is a positive value. The initial 

ar 

values of photon generation rate in Fig. 4 are positive. 
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FIG. 3. Photon generation rate for squeezed and coherent states for real value of a. 



FIG. 4. Photon generation rate versus the dimensionless time r for vacuum state and squeezed 
vacuum state with various squeezed parameters (|£| = 0.8, |£| = 1, |CI = 1-2 and |£| = 1.5 , = 0.) 

IV. THE PHOTON ABSORPTION RATE IN NONLINEAR MEDIA 

Gong in his paper [15] derived a Kerr-medium-like Hamiltonian in the vacuum induced 
by the radiation pressure for optomechanical systems which shows the underlying nonlinear 
mechanism more clearly than other approaches. In this mechanism the third-order suscep¬ 
tibility can be obtained as in usual Kerr media. 

So we can simulate cavity with moving mirror as a Kerr-medium-like with non-linear sus¬ 
ceptibility. With regard to the relationship between the cavity with moving mirror and 
nonlinearity it is needed that the effect of light on nonlinear media is explored. Therefore 
in the present and next sections we investigate the effects of squeezed and coherent light on 
a nonlinear media with third-order susceptibility. 

The multiphoton absorption rates are, first, studied in semi-classical approach [34], At 
this scheme, the electric field is considered classical and the atomic wave function is the 
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solution of the time-dependent Schrodinger equation, 


= ( 37 ) 

where 

H = H 0 + V, (38) 

here H is the Hamiltonian operator of the system. H 0 is the Hamiltonian for a free system 
and V is the interaction energy in the form, 

V = -dE(t ), (39) 

where d = er is the electric dipole moment, f and —e is the displacement operator and the 
charge of electron, respectively. To find the general solution for the system, the standard 
ansatz is made: 


^(r,t) = Yhi{t)ui{r)e iu,t , (40) 

i 

here hi(t) gives the probability amplitude and ui(r) satisfy the eigenvalue equation, when 
the interaction is not present, 

H 0 ui(r) = Eiui(r) E t = (41) 

The eigenstates constitute a complete orthonormal set satisfying the condition 

J ui(r)u* m (r)dr = 5i m . 

Therefore the Schrodinger equation is written in the following form, 

= Y hiMVrraWe-**^, (42) 

i 

where c Ui m = uii — LO m and V m i are the matrix elements of the interaction V, 

V ml = J u* m (r)V Ul (r)d 3 r. (43) 

Eq. (42) can be solved by using the perturbation theory. The Hamiltonian Eq. (38) is 
replaced by 

H = H 0 + XV, (44) 
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here A is expansion parameter. Also, the hi(t) in Eq.(42) and ij)(r,t) in Eq.(40) can be 
expanded in powers of A: 

hi(t) = + A+ A% (2) (f) + ..., 

ip(r,t) = -0 (o) (r, f)(f) + A ^ (1) (r,f)(f) + A 2 -0 (2) (r, t)(t) + ..., (45) 

here /i 7 ° corresponds to the ground state in the absence of the applied field and we set it 
equal to 1. So ip(°\r,t) can be considered as, 

ip(°\r,t) = u g (r)e~ iLJgt , l = g (46) 

By equating powers of A on each side in Eq. (42) the set of equations for different orders of 
hi(t) and h m (t) are obtained. 


dh m(t) 

dt 


~-(ih) 1 Y. h\ N l \t)V m i(r)e 


N = 1,2,3,... (47) 

For simplicity the electric field is considered as a monochromatic wave of the form: 


E(t) = Ee~ iut + c.c. (48) 

The N —photon absorption rate for the atom in state m at time t is given by 

C = < 49 > 

where P^\t) = \hm\t)\ 2 is A - —photon absorption probability. Now let us to consider a 
three level atom with ( g,m,n,l ) states( see Fig.5). It is straightforward to calculate the 
transition rate for third photon absorption in semi-classical approach[34]: 

d gm d mn d nl E 3 




l 2Ttpf(uj L9 — 3a;). 


(50) 


H 3 (u mg - u)(u ng -2u) 

Here Pf(oj) denotes the density of final state, as a distribution on the frequency range. d 9m 
is the electric dipole transition moment between ^(ground state) and state m. Also the 
polarization for this system is derived by applying the third order perturbation, 


(p( 3 )) =(-0(° ) |d|'0( 3 )) + ('0^ 1 )|(i|'0( 2 ))+ 

{'0 (2) |rf|'0 (1) > + <'0 (3) |^|'0 (O) >- 


( 51 ) 
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and by using Eqs. (40), (43),(45) and (46) one gets 

(p (3) ) =^hf\t)d 9 \r)e- i “ l ° t + 

i 

h { n\t)h*£\t)d mn (r)e~ iw ™ t + c.c. (52) 

nm 

For simplicity it is assumed that virtual levels (Z, m, n) are nearly coincident with real atomic 
levels. On the other hand, by using second quantization, the average of polarization oper¬ 
ator can be obtained. At this scheme, by introducing Dirac notation for up(r) as | f3), the 
orthogonality and completeness relation take the forms, 

m') = v 

XI W)V\ = !• (53) 

P 

The following identification is well understood, 

\p){0\ = ( 54 ) 

here operator c^'(c^) annihilates (creates) electron at the state \/3')(\/3)). So that we obtain 

the form of the polarization operator p = er, 

p=J2 \/ 3 )(/ 3 \ ef; \p')(p'\ = ( 55 ) 

pp' w' 

This form of polarization is general and includes all possible situations that may be occured, 
so we choose the part that is related to the third order of polarization, 


<p«>) = #‘{ct 6,) + + C.C. 


(56) 


Comparing the Eq. (56) and Eq. (52), we find that hf(t)e lU}3lt in semi-classical approach is 
equivalent to (c^q) in quantum approach, 

|fc< 3 >(t)| 2 =|<cjc,)| 2 . (57) 


Now, we need to evaluate the dynamics of microscopic polarization operator ( p l9 ) = 


[c]c g ). 


V. QUANTUM APPROACH 

The system under study is schematically shown in Fig. 5. 
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FIG. 5. Schematic representation of the third-order polarization. 

The total Hamiltonian for this system, includes two main parts, 

H tot = H 0 + H d , (58) 

here H 0 includes the non-interacting parts of carrier and photons system. The dipole inter¬ 
action between light and atom is given by H D . The explicit form of Hamiltonian in second 
quantization scheme as is given in [38-40], 


H 0 =£ a c\c g + £ m c f m c m + £ n c f n c n + £ l c\ci 
+ hw p (blbp + -) + fiWqftqbg + -) 

-\-ftjjj r (b\b r + -), 

h d = - {(ihf p d m9 clc g bp - ihr v d? m c\c m b\) 

+{ihT q d nm cic m b q - ihF q d mn $ m c n b\) 

+{ihT r d ln c\c n b r - ihF r d nl c\c x b +)}, (59) 

where T v is the interaction strength between carriers and photon with mode indx p. b p {b £) 
is the photon annihilation (creation) operator. By using the Heisenberg equation of motion 
the dynamics of the expectation value of polarization is given by 

d 

lh Q- t ^ 9 ) = ([c]c g ,H tot ]), 

ih ^(4 d a) = ( £9 ~ £l ){c\c g ) - { —ihT*d 9m (c\CjJJp) 

+ ihJ r *d nl (cl l Cgbl))}. (60) 
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In order to solve this quantum excitation, we have to evaluate the dynamics of the photon- 
assisted polarization. Applying the Heisenberg equation of motion for {c\c m b^ leades to: 


ih^(c\c m b ] p ) = (- e l + e m - hw p )(c\c m bl) 

+ ihF* q d nm *(c\c n blbl) 

- ihFpdT 9 (c\c g b p bl) 

-ihr;d ln \$ n c m b\bl) 

- iUT v d m9 (cjcrnC^Cg ). (61) 

and for ($ n c g b]) one obtains: 

ih ^l&g b l) = (-£ n + £ 9 - hu r )(clc g bl) 

+ ihF;d m °*(cl l cJtbl) 

+ ihF r d ln {c\c g b r b\) 

- ihr q d nm \c*j. g b\bl) 

+ i?iJ : r d ln {$ l c n $ n Cg). (62) 


Eq. (61) and Eq. (62) have the same intuitive interpretation. In Eq. (61), the photon- 
assisted polarization is generated either via spontaneous emission or via absorbing a new 
photon, including the (cjCnbfy) and ( c\c g b p b j,) terms, respectively. It may be also transform 
to a new term, via spontaneous emission by (c^c m 6 j, 6 j) term. The latter term contains the 
correlated carriers. For our purpose, it is sufficient to consider only the terms which produce 
photon-assisted by spontaneous emission. So we consider the other terms as decay terms. 
The dynamics of these terms follow from Heisenberg equation of motion, 

ih^{c\c n blb\) =(-e l + e n - hw p - hw q ){c\c n blb\) 

-ihF r d ln *r(blblb\) 

+ ihF r d ln * f l {btblb\) 

-i^i(c\c n b ] p b\), (63) 

where f n is the density of state, n and 711 act as dephasing source with elements defined in 
appendix A. The other spontaneous terms, (c\c m b\w p ) and (c^CgwJ)]) in Eqs. (61) and (62), 
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are structurally similar to Eq. (63) 


d 


ih— (c ] n c m b\bl) ={-e n + e m - frw p - hw r )(clc m blbl) 

- ihr q d nm * f m (b\b\bl) 

+ ihF q d nm * f n {b\b\bl) 


- i^\2{c' n c m blb ] p ). 


t, 


( 64 ) 


and 


d 


ih di^ £ gW =(~ £m + e 9 - hw q - hu r )^ m c g b\bl) 

+ ihr p d m9 * r(blb\bl) 

- ihr q d m9 * f 9 (b\bibi) 




- t^22{C' m C g 0' q 


$$)■ 


(65) 


The explicit forms of 712 and 722 are given in appendix A. 


VI. ANALYTIC SOLUTION 

The expectation value for any operator 0 can be written according to the density matrix/?: 

( 6 ) = Tr(pO). ( 66 ) 

For our subsequent calculations it is advantageous to make the unitary transformation p = 
f/ _ 1 pf/, in it, 

U(t) = e - iH ° t/h , (67) 

where H 0 has been given in Eq.(59). Substituting p by p in Eq. ( 66 ), the expectation value 
yields, 

Tr(pO) = Tr{U~pU~ 1 0) = Tr(e- iHot/h pe iH ° t/h 6). ( 68 ) 

Considering the general property of the trace Tr(ABC) = Tr(CAB), 

Tr(e- iHot/h pe iHot/n 6) = Tr(pe iHot/h 6e- iH ° t/h ). (69) 
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Accordingly, for the operators appeared in the Eqs. (60)-(65), by using the following relation, 


a 


e~ aA Oe aA = O - a[A, O] + —[A, [A, O ]] + 


the expectation values can be calculated, 




b\b\) = V clc n O' p O' q )e 


i (e l — £ n -\-ftWp -\~hbJq ) t /fl 




(70) 


where 


(c\c n b ] p b\) = (, o(c\cJ>lb\)) 

(blblbl) = <p(S+SjS+)). (71) 


By inserting Eqs. (70) in Eq. (63), and considering the decay term as a exponential term, 
then differentiation with respect to the t, the solution will be given by: 



n u p u q 


- f n )(blb\bl) 

+ e n + huj r — ih'y n)) 

h 

[exp(^(—e* + £ n + hw r - ih'-fn)t) - 1]. 
n 


(72) 


And so on for Eqs. (64) and (65) we have, 


(Cn,C rn brbp') 


{CmCgbqbr ) 


ihJ r *d nm * (f n - f r 


ih(-(-£ n + £ m + hw q - ih^ i 2 )) 
n 

[exp (^(-£ n -j- £ m + fVjJq i^ 7l2 )t) - 1], 

ihF;dr°*{f m - f 9 )(blb\bt) 

ih(^(-£ m + £ 9 + hjjp - ih 722 )) 


(73) 


[exp(-(—+ e 9 + - ih^ 2 2 )t) - 1]. 


(74) 
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If we transform the expectation value in Eqs. (61) and (62) according to the Eq.(69) and 
inserting Eqs. (72)-(78) into these equations, 


(i cjc m b f p ) = 

ihF*d nm *ihF;d ln * (, f l - f n ) (blb\bt) 

(~£ l + £ n + huj r — *^7ll )(—£ l + £ m + hjj r + fto q — ih'yi) 

ihF*d nm *iHF*d ln * (, f n - f m ) (b ] p b\bl) 

(—£ n + £ m + hu q — + £ m + ftO r + ftO q — ihjl) 

(exp( —(—£ l + £ m + fto r + ftjq — ihrfi)t) — 1). (75) 

n 

and 

(cUgbt) = 

iHF* d nm *iHF*d m9 * ( f n - f m ) (SfSJSJ) 

(-£ n + £ m + ftjq - ih^i 2 )(-£ n + £ 9 + ftJp + hw q - ihrj-i) 

ihF*d nm *itu^cr 9 * (f m - f 9 ) (SJSJSJ) 

(-£ m + £ 9 + ftJp - ih^22)(-£ n + £ 9 + fto p + ftjq - itvyz) 

(exp(^(-£ n + £ 9 + ttjp + frw q - ih^ 2 )t) - 1). (76) 

By replacing these results into Eq. (60), the final form can be obtained. Now, we consider 
some assumptions for description the three photon absorption process. The states n and 
m are the virtual states, so their populations are vanishes. In the three photon absorption, 
in Fig 5, we assume the involved frequencies are the same, so that, co p = oj q = ui r = uj so 
F p = F q = F r = F and F q = F p = b\ = F. On the other hand the virtuale states are nearly 
coincident with one of the resonant real level. Therefore we consider fto ln = ftj nrn = ft:™ 9 . 
For simplicity the decay rate of states has been considered as the same to each other as 7 . 
According to this assumption, 

{c\c g ) = 

ihF* r d ln * ihF*d nm *iKF*pd m9 *if 1 - f 9 )(blb\bt) 
h 3 (u ln — uj + i^)(uj lrn — 2uj + i'j) 
exp (-i{u la - 3w + * 7 )*) - 1 

(uj 19 — 3u; + i'y) ’ 1 ’ 

According to Eqs. (69) and (70) one obtains, 

( 6 + 6 + 6 +) = (blb\bl)e i ^ r+nwp+hw ^ t/h (78) 
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Inserting the above equations into the Eq. (77), the suitable form for expectation value will 
be achieved, 

{c\c g ) = 

ihF;d ln *ihF*d nm *ihF;d m9 \f l - f 9 )(blb\bl) 

H 3 (u ln — u + i^){u> lrn — 2 u + i'j) 
exp (—z(o/ 9 — 3u + i^)t) — 1 
(a : lg — 3u + ij) 

^i{£ l —£ 9 )t/h^—i{hbj r +tiwp+hbj q )tlh ( 79 ) 

here hw lm = e l — e m . Using Eq. (49) and Eq. (57), one obtains the rate of three photon 

absorption in quantum approach, 

r3 (ihF*d ln *)(ihT*d nm *)(ihT*d m 9 *)(f l - p)(bW) 2 
1 h 3 (u ln — u + — 2 oj + 7y) 

x 2vrp(cn' 9 -3cn), (80) 


where p(uj l9 — 3o;) is the density of final state and Lorentzian shape is a well-known example 
of it and its derivation is given in Appendix B. In quantum approach the electric held has 
the form as E = E~ + E + and the negative part defined as E~ = ihE*b f . Therefore the form 
of the absorption rate in quantum and semi-classical approaches is the same (see the Eq. 
(50)). By considering Eq. (21) one can obtain the expectation value (b^b^b^) for coherent 
and squeezed state. The normalized rate of the third order absorption for squeezed and 


coherent lights with the same intensity, arc illustrated in Fig. 6, where 7 = 27t(1 x 10 13 ) 


Rad 
Sec ’ 

d ln = d nm = d mg = 8 x 10~ 30 cm, |ck | 2 = 7 and different values of the squeezed parameter 
have been chosen. 

This figure shows that nonlinear media are very sensitive to the state of light. When the 
incident light is in a quantum state, specially squeezed one, the rate of absorption is several 
times greater than the case if the coherent light is applied. On the other hand, one can see 
obviously that for the same intensity, by increasing the squeezed parameter the difference 
of the rate absorption between coherent and squeezed light becomes very prominent. 

This important point is an indication for the accepted result, that the nonlinearity can 
be used for amplification of non-classicality. So that the proposed system can be considered 
as a detector for the squeezed state of radiation. 

Investigating figures 1 and 6 shows that the photon generation by a driving held in 
squeezed state, in ID cavity with moving wall and three photon absorption in nonlinear 
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FIG. 6. The rate of three photon absorption for squeezed and coherent light with the same intensity 
and various squeezed parameters (|£| = 0.8, |£| = 1, |£| = 1.2 and |£| = 1.5 , = 0. Rm ax is 

related to the maximum value of squeezed light with |£| = 1.5 ). 

medium, is several times greater than the case of coherent light. It seems that this phe¬ 
nomenon may be explained by considering that the noises in squeezed light, at least in one 
of its quadrature, is suppressed, and altogether the squeezed light shows more interaction 
power in nonlinear media. 


VII. CONCLUSION 

The main purposes of the paper are twofold. First, we have studied one-dimensional cav¬ 
ity with moving mirror in presence of squeezed and coherent light. The photon creation rate 
has shown an obvious difference between these two modes of light. The photon generation 
rate increased with increasing of squeezing parameter. 

Second, by another approach, we have demonstrated that the nonlinear responses of the 
media is highly affected by the state of the interaction field. We have illustrated the third 
order of the photon absorption rate of squeezed light in comparison to the coherent case, 
dependent on the squeezing parameter, indicates drastic difference. 

As most important result, we conclude that the nonlinear effects in ID cavity with moving 
mirror are comparable with a nonlinear medium, when they interact with different states of 
light. In other words, the nonlinearity amplifies the non-classicality of states. We have seen 
the drastic difference between squeezed and coherent states of light, which can be used for 
detection the non-classicality of states. 
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Appendix A 


The decaying terms in Eqs. (72)-(76) have the following forms, 

h n = - ihF q d nm (c\c m b q blb\) 

- ihF p d mg (c ] m c g c\c n b\) 

- ihF q d nm ($ n c m c\c n bl >, (Al) 


*712 = - ihT r d ln {c\c m b r b ] p bl) 

+ ihF p d m9 (cl c g b p blbl) 

- ihl'rd 171 (c]c n clc m bl) 

- ihFpd™ 9 (c ] m c g c ] n c m bl ), 


and 


?722 = - ihT q d nm {$ n c g b q b ] q b\) 
-ihT q d™{cic m clcgb\) 
- ihJ r r d ln (c]CnclnCgbg). 


(A2) 


(A3) 


Appendix B 

According to Ref. [34], the time dependence of the square of absolute value of the 
expectation value can be obtained, 

exp(— i(uj l9 — 3co + ry)t) — 1 2 4 sin 2 ((o / 9 — 3a; + i^)t/2) 

(oj lg — 3u + i'f) (■ oj 19 — 3u + i'y) 2 

On the other hand, for large values of the time, one can prove, 

lim^noi 2 ——= 2irt8UJ 9 — 3oj) for x = ( u l9 — 3u)t/2. (B2) 

x l 

In real situation, the final state is not sharp in the frequency transition and make a contin¬ 
uous distribution in frequency domain. So the 5(u> l9 — 3u) can be replaced by p[oj l9 — 3a;). 
Usually, the p{pj l9 — 3 uj) has Lorentzian shape including the decay term 7 . So the absorption 
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probability can be written in the form, 

p3 ihF*d ln *ihF*d nm *ihF*d m9 \f l - f 9 )(titi&) 2 

1 h 3 (oj ln — to + — 2oj + i'y) 

x 27 Ttp{uj l9 — 3a;). 

Using Eq.(49) one gets the Eq.(80). 
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